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Abstract 

Isolated domains in lipid monolayers exhibit shape transitions as 
a consequence of the competing effects of line tension and long-range 
dipolar interactions. Nontrivial integrals ultimately arise in stability 
analysis calculations of these domains. In this work we present a closed 
form evaluation of such integrals in terms of associated Legendre func- 
tions. Our approach is motivated by and extends the work of J. M. 
Deutch and F. E. Low (J. Phys. Chem. 1992, 96, 7097-7101). Our 
closed form solutions lead to simpler analytic expressions for stability 
thresholds, which are easier to evaluate explicitly. Relevant asymp- 
totic behavior is calculated in various limits of interest. For the case 
of Coulomb interactions, a general closed form expression is obtained 
for the critical domain radii. 
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1 Introduction 

In the past few years the study of shape transitions in hpid monolayer do- 
mains at the air-water interface has received considerable attention from the 
physical-chemical community The great interest towards this subject 
is in part due to the diversity of domain shapes, which range from simple 
circular structures to complex labyrinthine patterns. Theoretical investiga- 
tions study this variety of shapes by a model in which an array of molecular 
dipoles is arranged normal to the two-dimensional domain. The basic mech- 
anism behind the formation of such structures arises from the competition 
between interfacial tension and long-range electrostatic interactions. Line 
tension acts to minimize the domain perimeter and favors compact domain 
shapes. In contrast, repulsive electrostatic dipole-dipole interaction tends to 
maximize the distance between dipoles and favors extented domain shapes. 
In this context, the total area of the domain is not changed. 

In previous work McConnell p[, Deutch and Low 0, and Goldstein and 
Jackson 0] have studied the shape transitions of initially circular (radius R) 
lipid monolayer domains in a lowest order (quadratic) stability analysis. In all 
these theoretical studies many relevant quantities are calculated exactly, but 
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due to their complexity, some of the computed expressions are not evaluated 
in closed form. For example, the formulas for stability thresholds found 
in 1^, D, Q are usually expressed in terms of nontrivial integrals. Related 
integrals arise in stability investigations of other systems, also presenting 
polarized domains, such as planar- confined ferrofluids in magnetic fields 
and thin ferromagnetic films forming magnetic bubble patterns 0. 

A considerable difficulty to perform stability analysis of lipid monolayer 
domains is the evaluation of the electrostatic dipole-dipole energies for var- 
ious domain shapes. These calculations involve the manipulation of double 
line integrals which are obtained by an ingenious application of Green's the- 
orem [0. The subtleties concerning the validity of Green's theorem under 
such circumstances are connected to the necessity of introducing a cutoff dis- 
tance required to prevent divergence at zero separation between dipoles [§]. 
As shown by Deutch and Low [Q], the subtle application of Green's theorem 
can be avoided by calculating the electrostatic energy directly through area 
integrals. 

The purpose of the present work is to evaluate in closed form the integrals 
and pertinent stability conditions obtained in reference 0|. Even though 
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the nontrivial integrals which appear in references [0, Q basically belong 
to the same general family, we focus on Deutch and Low's paper since 
their work precisely illustrates the utility of our closed form expressions. 
Our calculations are performed in the same spirit of those by Miranda and 
Widom 0, who recently carried out a related closed form study focused 
on ferrofluids. We evaluate the integrals in terms of associated Legendre 
functions. Asymptotic behavior in interesting limits such as small cutoff, or 
high spatial frequency are rigorously calculated. Finally, we address the case 
which considers repulsive Coulombic interactions within a charged domain. 
Comparison between our results and those derived in is consistently made 
throughout the text. We obtain simpler and compact expressions which 
facilitate the explicit evaluation of many important quantities derived in 
previous studies of shape transitions in lipid monolayers P, H 

2 Stability analysis of lipid monolayer domains 
- Closed form results 
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2.1 The integrals Bn{a) 

Following Deutch and Low the perturbed domain shape is described by 
the Fourier expansion 

r(^) = i?+^CnexpM), (1) 

where 6 is the polar angle, n (an integer) denotes the discrete azimuthal 
wave number and (n represents the perturbation amplitudes (figure 1). In 
equation ([1|) (q is expressed in terms of the amplitudes Cn ^ 0) in order 
to satisfy the condition that the total A = (1/2) f^^'r {6 fde of the 
domain remains constant and independent of the perturbation, i.e., A = tcR"^. 
The perturbation mode with \n\ = 1 corresponds to a global off-center shift 
(rigid translation) of the unperturbed initial shape. It is a spurious, shape- 
preserving mode that does not contribute to domain energy changes. 

The total energy of the domain described by r{9) can be calculated by 
considering electrostatic and line-tension contributions. The bulk (dipole- 
dipole) electrostatic energy of a distribution of molecular dipoles oriented 
in a direction normal to the two-dimensional domain is given by the double 
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surface integral 
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where is the dipole density of the surface and a is a cutoff distance of 
closest approach of neighboring dipoles. On the other hand, the line-tension 
energy, given by the product of line tension A by the domain perimeter, is 



Fx[r{e)] = x r de 

Jo 



r{9f + 
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\ de ) 
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By using equations Deutch and Low showed that the total en- 

ergy F[r{6)] = W[r{9)] +Fx[r{9)], accurate to order C(C^), can be expressed 



in the form 

^ n>0 

where Fq is the energy of the unperturbed circle and 



(4) 



= X{n^ - 1) - 2/i2[5i(d) - 
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where 



cos n6 



3/2 



de 
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"(l + f) -cos^ 

and a = a/R. We can examine the stability of the circle simply by looking 



at the sign of Qn'- instability occurs when f2„ < 0. Similar results have been 



obtained in references and |Q, despite subtleties regarding their implemen- 
tation of the cutoff. Even though expression (H) is exact (to second order), 
it is written in terms of the integrals and have not been explicitly 

evaluated. Consequently, relevant physical quantities like the critical radii 
Rn (radii beyond which a circular domain becomes unstable with respect 
to n-fold harmonic perturbations) are often written in an indirect fashion, 
which involves complicated sums and integrals, or requires the use of recur- 
sion relations [|l^ . The lack of closed form results makes difficult the precise 



and explicit evaluation of Rn for arbitrary n. Therefore, it is of practical 
interest to find closed form solutions for the integrals i?„(d;). 

The closed form expression for the integrals Bn{&) can be found with the 
help of the integral representation for the associated Legendre function of 
the second kind [jll| 



COS z/vr 



oo 



exp 



[2 + cosht]™+i/s 



dt\, (7) 



with Re(m) > —1/2, Re(z/ + m) > — 1 and |arg(z ± 1)| < tt. T denotes the 
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gamma function. Set v = n — 1/2 and m = 1 in equation (|^) and compare 



the resulting expression with equation to get 




QU/2 (i + f ) 




(8) 



This way, equation (|^) can be rewritten as 



rtr, = A(ra2 - 1) - 




(9) 



We stress that n„ is now written in a closed form fashion, which explicitly 
expresses its precise functional dependence on n, a, X and fi. The fact that fin 
is written in terms of the associated Legendre functions is very convenient, 
since their fundamental properties and basic functional behavior are quite 



in the calculation of asymptotic expressions in important limiting situations 
such as small cutoff (a — > 0) and long wavelengths (n ^ oo, a — > 0, with 
a kept finite and small). In such relevant limits our closed solution leads 
to simple and elegant expressions, which are easier to evaluate than those 
derived in 0. 

In order to gain insight about the role of fi„ in the domain instability, 
we plot in figure 2 the ratio fi„/2A as a function of mode number n. We 



well known ||TT|1. The usefulness 




specially evident 
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assume that d = 2 x 10^^ and take three different values of the dimension- 
less parameter /x^/A, a quantity that characterizes the relative importance of 
dipolar and surface energies. All points are obtained exactly, using our equa- 
tion (P). The domain boundary remains stable to n-th mode distortions up 
to a critical value Dcr{n) of /x^/A, defined by setting fi„ = in equation (P). 
In (a) yU^/A < Dcri!^), so i7„ > and the domain is stable. For larger values 
of /x^/A (see (b) and (c)) f2„ can be negative and the domain deforms. 

2.2 Stability condition and critical radii Rn 

For amphiphilic systems the ratio between the microscopic cutoff a and the 
radius R of the circular domain is a small quantity. In this section, we take 
advantage of our closed form expression (|) to evaluate the stability condition 
and critical radii Rn in the limit of small a. 

From our discussion in the previous section the stability condition i7„ > 
for circular lipid monolayer domains can be written as 




2 (i + f ) - Q\/2 (i + f ) 



(10) 



«(i + f)' 



In order to evaluate (|T0|) in the limit of small a, we use the following property 
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of the associated Legendre functions [|n| 



(11) 



to rewrite (|TOD as 



g„_i/2(i + 2x) -gi/2(i + 2x) 



(n2 - 1) 



(12) 



where x = 



To get a closed form expression for the stabihty condition (|T2[) we need 
to evaluate (5n-i/2(l + 2a;) for small values of x. To succeed in doing this, 
we use reference |jl2[ to expand Qn- 1/2(1 + 2x) near x = 0, to first order in x 



Qn-i/2il + 2x) 



-C-llnf-^ 

2 \x+l 



ilj{n+- 



1-C) 
+ 0{x^ lnx,x^), 



2 \x + 



X 



(13) 



where Euler's psi function is the logarithmic derivative of the gamma 



function |Tl| and C is Euler's constant. 

Using equations (p!^ ) and (p!3|), we obtain the closed form stability condi- 



tion, accurate to order a, 
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At this point, we compare our expression (|I^) with an equivalent result ob- 
tained by Deutch and Low p| (equation (36) in their work). In reference 
the part corresponding to the term between curly brackets on the right hand 
side of ( |14|) is written as complicated definite integrals involving Chebyshev 
polynomials. The authors in do not solve such integrals explicitly. In con- 
trast, our derivation, based on the well known properties of the associated 
Legendre functions, naturally provides a closed form stability condition. Ob- 
viously, our solutions are useful only if they lead to easier manipulation of the 
obtained results. This is exactly the case: the right hand side of our stability 
condition (|1^) is much easier to be evaluated explicitly than the equivalent 
piece which appears in Q]. Notice that analytic evaluation of equation (0) 
can be promptly performed, for specific n, by means of a simple property 



relating ip functions ITT 



H"4)-K"4) = 2S^^ (15) 

Another important quantity that can be immediatly obtained from (0) 
refers to the critical radii at which the circular domain becomes unstable 
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with respect to a given Fourier mode n 



i?n = ^ exp ( ^ ) exp (Z„) , (16) 



where the critical paramenter for transition is given by 



_ 3 1 Un^ - r 



2/ " V2 



(17) 



The values calculated for Z„ directly using our equation (0) agree precisely 
with those computed in equation (39) of reference [Q. Again, the major 
difference between our results resides on the fact that in to get a particular 
value of Zn it is necessary to appeal to an indirect approach which requires 
solving several different integrals involving Chebyshev polynomials, for each 
n to be considered. Our equation ( [T^ ) offers a simpler, direct and more 
convenient way to evaluate the various Zn- 

Figure 3 depicts the dimensionless critical radius Rn/o. plotted against 
the ratio A//i^. All curves are computed analytically from our equations (p^) 
and (p!7|). This figure illustrates important physical behavior captured by 
lowest order stability analysis. Mode n = 2 is always the first to become 
unstable (elliptical instability), establishing a stability boundary. At fixed 
value of A//i^, higher modes (n > 3) become unstable at larger Rn/a- As 
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Rn/ot increases, each mode becomes dominant for a range of radii, and hence 
there is a cascade to higher modes. For a given critical radius, smaller values 
of A//i^ lead to stronger mode competition, meaning more intrincated domain 
boundaries. Related stability diagrams have been numerically obtained in P| 



and ||T3[. The critical parameter Zn is plotted as a function of mode number 



n in the inset of figure 3. 

2.3 Relevant asymptotic behavior and the Coulomb in- 
teraction case 

We begin this section calculating further asymptotic expansions by taking 
advantage of our closed form expression (|TI1|). First, consider the large n and 
small a limit, assuming that the product na is kept fixed. This is the limit of 
interest in many practical situations p|, |^, ^. Considering the large n limit. 



use the asymptotic expression 



[cosh (^) 



Km{r), (18) 



n^oo ?T,'"exp(^?7rmj 
where Km{T) is the modified Bessel function of order m, and set 



T = n cosh ^ (l + ^) ^na + 0{na^), (19) 



2 
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to rewrite stability condition (|T0|) as 



A ^ 2 



Ki{a) — nKi{na) 



a 1 + 



1/2 



(20) 



In addition, take into account the fact that a is small, keeping na fixed, to 



obtain 



> 2 



mar na 



-Ki(na) 



(21) 



where we have used the small argument expansion 



1 1 
-y^2 



y + 0{y^ lny,y^ 



(22) 



and a power series expansion of the denominator in (pOl) . Equation (pTD is 



quite general, being valid for large values of n and small values of a, provided 
na is kept finite. We would like to point out that our result (pT]) differs from 
the equivalent one obtained in ||^ (see equation (46) in their work). Because 
of a mistake in their calculation, Deutch and Low have written their equation 
(46) in terms of the modified Bessel function of order zero Ko{na) (there is 
also a sign error), which is not the correct result. This minor mistake would 
be of no consequence if the authors in |^ would have not derived an important 
subsequent result, using their incorrect expression (46) as a starting point. 
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This fact generates some confusion. Further comments on this issue are given 
in the next paragraph. 

We take advantage of our equation (^l]) to discuss the so-called straight 



edge limit p, [Ij]. This is the situation in which the domain interface is 
flat with a small perturbation of wave number k. To get this limit from our 
circular geometry, we take radius i? ^ oo at fixed a. Fixed wave number k = 
n/R implies n —>■ oo and a —>■ 0, with fixed na as in our asymptotic expansion 
equation (0). Taking the limit of (|21|) for large wavelengths compared to 
the cutoff {ka <C 1), we use the small argument expansion (^) to get the 
critical spatial frequency 



a fx ^ 1^' ' 
- exp I — + C - - 



(23) 



2 ' 2^ 
The straight edge of the domain is unstable for all spatial frequencies k < kc- 



This result agrees with equation (13) in McConnell's work [|T4| since the line 



tension considered in his derivation contains an additional electrostatic con- 



tribution given by — /i [|T5|. Our result p^) also agrees with the quantity 
i?7v/(n) = n/kc derived by Deutch and Low (equation (48)). In principle, we 
should not expect agreement since RM{n) in is obtained using their incor- 
rect intermediate expressions (46) and (47), as pointed out in the previous 
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paragraph. It is worth mentioning that the critical spatial frequency (p^) 
could have been calculated by taking the large n limit of our earlier re- 
sults (|1|) and (|17D- 

We conclude by extending our analysis to the case of Coulomb interactions 



in charged domains H^]. As indicated in 0, for a Coulomb interaction 
between two particles (charge q) of the form u{p) = [q'^/p), with p = [r^ + 
_ 2rr'cos^^]^/^, the domain stability condition can be written as 

X{n' - 1) > 2'/'q'R'[B^ie) - B^e)], (24) 

where 

B„(,) , r '"''"l de. (25) 
•^0 [1 — cosb'J ' 

As in the case of the dipolar interaction potential, the integrals Bn{9) are not 
evaluated in 0. Using equation (|^) with u = n — 1/2, m = Q and z ^ 1, plus 
expansion (0), we find closed form expressions for the integrals Bn{9), and 
consequently for the stability condition (p^. With the help of these results, 
we calculate the critical radii for a given n mode distortion 



Rr(n) 



1/2 



2g2 [V^ (n + i)-^(f) 



(26) 
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Keller et al. [jl6[ studied the specific case of elliptical deformations in charged 
domains. Our result is identical to theirs if we set n = 2 in equation (|26|). 
Equation (^) generalizes previous results found for Rc{n) [§, providing 
a closed form expression, valid for arbitrary n. 

3 Concluding Remarks 

Lipid monolayers form two-dimensional domains as a result of the compe- 
tition between line-tension and long-range repulsive interactions. Stability 
analysis of these domains depends on a family of nontrivial integrals. For 
repulsive electrostatic dipole-dipole interaction, we present a closed form eval- 
uation of these integrals as a combination of associated Legendre functions. 
Asymptotic behavior in several interesting limits are evaluated. A general 
closed form result is also obtained for the case in which the repulsive inter- 
action is of Coulombic nature. Our theoretical results should be viewed as 
complementary to the work of Deutch and Low leading to exact, sim- 
pler and explicit formulae which commonly appear in the study of shape 
transitions in lipid monolayer domains. 

We hope our results will be helpful for future work on the subject of 
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shape instabilities in amphiphilic monolayers. While the present work has 
addressed transitions associated to a slightly deformed circle, we anticipate 
that the analytical treatment of more complicated domain shapes may bene- 
fit from our closed form results. In higher order perturbative calculations, the 
various harmonic modes couple and cannot be studied separately. In order 
to cope with larger distortions a mode coupling theory is required. We point 
out that a weakly nonlinear analysis of such situation can be carried out ana- 
lytically, following a related recent work on viscous fingering instability [|17 



Once applied to monolayers, this approach would also lead to nontrivial in- 
tegrals which could be explicitly evaluated by employing the methodology 
and results we presented here. 
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Figure Captions 



Figure 1: Schematic configuration of a two-dimensional lipid domain. The 
dashed line represents an initially circular domain of radius R and the solid 
undulated curve depicts its perturbed shape r{0)= R + (■ The perturbation 
C= J2 CnGxp {in9), where n denotes the azimuthal wave number and 9 gives 
the angular location of the points on the boundary. The domain area enclosed 
by the curve r{6) remains constant and equals to ttR"^. 

Figure 2: Ratio r2n/2A as a function of mode number n, for three distinct 
values of fi'^/X (a) 0.090; (b) 0.107; (c) 0.112. The critical dimensionless 
parameter Dcr{'^) ~ 0.095 and a = 2 x 10~^. All points were found exactly 
using our closed form equation (|^). 

Figure 3: Dimensionless critical radius Rn/a as a function of A//i^. Closed 
form equations ([TBI) and (^) were used to plot the various curves 2 < n < 16. 
Inset: behavior of the critical parameter Z„ with respect to n. 
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